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We present a study of the experimental determination of the forward-backward asymmetry in the 
process e+e" -s- tt and in the subsequent t -)■ Wb decay, studied in the context of the International 
Linear Collider. This process probes the elementary couplings of the top quark to the photon, the Z 
and the W bosons at a level of precision that is difficult to achieve at hadron colliders. Measurement 
of the forward-backward asymmetry requires excellent b quark identification and determination of the 
quark charge. The study reported here is performed in the most challenging all-hadronic channel 
e+e" -^ bbqqqq. It includes realistic details of the experimental environment, a full IVIonte Carlo 
simulation of the detector, based on the Silicon Detector concept, and realistic event reconstruction. 
The forward-backward asymmetries are determined to a precision of approximately 1% for each of 
two choices of beam polarization. We analyze the implications for the determination of the tiZ and 
Wtb couplings. 

PACS numbers: 14.65.Ha 



I. INTRODUCTION 



Tine top quark is substantially more massive tinan tine 
otiner known quarks. Simply by virtue of tinis fact, tine top 
quark couples more strongly to the particles tinat gen- 
erate tine spontaneous symmetry breaking of tine elec- 
troweak interactions. It is possible tinat tine large mass of 
tine top quark is explained by new interactions of tine top 
quark. It is tinus important to measure tine properties of 
tine top quark carefully, searcining for signals of special 
interactions of this quark. 

Particularly interesting quantities to study are the form 
factors that describe the coupling of the top quark to 
elementary currents. These are the analogues for any 
new interactions of the proton form factors, which played 
such a large role in the elucidation of QCD. We will study 
the process e+e^ -^ ti. In this reaction, two sets of form 
factors enter, the form factors that describe the 7 and Z 
couplings to ti, which describe the tt production vertex, 
andthe the form factors that describe the W coupling 
to th, which describe the t and t decay vertices. As a 
matter of principle, a full reconstruction of the tt sys- 
tem in e+e" annihilation can give information on both 
sets of vertices. The effects of the possible form fac- 
tors on observables of the tt system have been studied 
by many authors, for example,^"''. Some of these cou- 
plings will be constrained at the Large Hadron Collider 
(LHC), but others are very difficult to access there. In 
particular, the vector and axial vector couplings of the 
top quark to the Z boson are shifted by new physics ef- 
fects in many models^i"— . However, these couplings are 



very difficult to measure precisely at the LHC, and the 
associated form factors are completely inaccessible at 
values of Q^ larger than m|. 

In this paper, we will begin a study of the determina- 
tion of these form factors under realistic experimental 
conditions at the proposed future e+e^ collider, the In- 
ternational Linear Collider (ILC). We will make use of 
the detector model given by the Silicon Detector (SiD) 
concept and the set of full-simulation tools developed for 
the benchmarking of SiDH. These tools provide a very 
detailed simulation of the experimental environment at 
the ILC. 

We will consider the forward-backward asymmetries 
both for the 5 and b quarks and for the t and t quarks. In 
each case, the forward-backward asymmetry is defined 
as: 



Afb = 



a{e < 90°) - cr{e > 90°) 
cr(6l < 90°) + cr(6' > 90°) 



(1) 
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where a{6 < 90°) is the cross section of the events in 
which the bor t quark has a polar angle of less than 90° 
in the centre of mass frame of reference. The standard 
spherical coordinate system convention is used to define 
9. This asymmetry measurement is a complex analysis 
in a dense multi-jet environment. Typical events have 6- 
jet final states. Flavor-tagging must be done to identify 
the b quark jets and resolve the combinatoric ambigu- 
ities. Quark charges must be measured to distinguish 
the t and i decay products. Detector resolution and 
acceptance together with non-ideal efficiency and pu- 
rity of the reconstruction algorithms could play a critical 
role in determining the ultimate sensitivity of the mea- 
surement and hence its physics reach. This study ad- 
dresses these issues for the first time. Our conclusion is 
that, with the beam conditions and integrated luminosi- 
ties that the ILC will provide, a well-designed detector 



can overcome these potential problems and realize the 
small measurement uncertainties that were projected in 
parametric studies. 

The paper is organized as follows: Section II gives 
general parameters of top quark production at the ILC. 
Section III introduces the SID detector concept. Section 
IV presents the software framework used in this anal- 
ysis. Section V describes the signal selection and the 
calculation of the cross section for the fully hadronic tt fi- 
nal state. This section also discusses the flavor-tagging 
method and its performance. Section VI is devoted to 
the quark charge reconstruction algorithms, which are 
fundamental to the analysis. The results for forward- 
backward asymmetries are presented in Section VII. 
Section VIM puts these results in context by interpret- 
ing them as bounds on deviations of the ZU and Wtb 
form factors from their Standard Model values. Section 
IX gives out conclusions. 
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FIG. 1 : Disposition of subdetectors in SiD quadrant. All dimen- 
sions are in mm. 



II. TOP QUARK AT THE ILC 

The International Linear Collider is a proposed 
electron-positron accelerator operating in the centre of 
mass energy range ^/s = 200 GeV - 500 GeV An up- 
grade to the centre of mass energy of 1 TeV is also en- 
visaged as are possible calibration runs at the Z boson 
mass energy"'^. The maximum design luminosity is 2 x 
10^^^ cm^^ s^^ In the analysis presented here the cen- 
tre of mass energy and total integrated luminosity were 
assumed to be respectively 500 GeV and 500 fb"\ the 
latter one equivalent to a few years of ILC running. 

The top quark at the ILC, assuming the 500 GeV op- 
eration, is mainly produced in pairs through the e+er -^ 



Z ^ tt and e+f 



7 



tt processes. The theoreti- 



cal total cross-section of top quark pair production is ap- 
proximately 600 fb^^. Although this value is substantially 
lower than the one at the Large Hadron Collider (LHC), 
the clean environment, well defined initial state and po- 
larization make the ILC an ideal machine to perform top 
quark precision measurements. 



III. THE SID DETECTOR CONCEPT 

The top quark properties are studied with Silicon De- 
tector concept which is a general purpose detector de- 
signed to perform precision measurements and at the 
same time to be sensitive to a wide range of possible 
new phenomena^ \ It is based on a five layer silicon pixel 
vertex detector, silicon tracking with single bunch time 
stamping capabilities, silicon-tungsten electromagnetic 
calorimetry and a highly segmented hadronic calorime- 
ter. The Particle Flow Algorithm (PFA)"''* is an important 
strategy driving the basic philosophy and layout of the 
detector. SiD also incorporates a five Tesia solenoid, 
an iron flux return and a muon identification system. A 
schematic view of SiD quadrant is shown in FigurelU 



IV. ANALYSIS FRAMEWORK 

The event generation has been performed using the 
WHIZARD MC generatori^ii^. Event samples were cre- 
ated with the expected ILC baseline parameters of 80% 
electron and 30% positron polarization. Half of the event 
sample was created with a positive electron and nega- 
tive positron polarization, while the other half has been 
created with a negative electron and positive positron 
polarization. 

WHIZARD was used to generate samples of all 0, 
2, 4, and 6 fermion final states as well as top quark- 
dominated 8 fermion processes. This generation used 
electroweak vertices only, with gluon emission turned 
off. The intent of this strategy was to correctly de- 
scribe multifermion processes such as return to the Z 
(e+e" ^ 7*Z* ^ 4 fermions) and similar processes with 
intermediate off-shell W bosons. These reactions with 
t-channel exchange and off-shell electroweak bosons 
are the most important backgrounds to multi-fermion 
e+e~ annihilation processes. QCD was included in the 
events by using PYTHIA to evolve final-state quarks 
through parton showering, fragmentation, and decay. 
PYTHIAJI was also used to generate final-state photon 
radiation. There is no double-counting of multi-fermion 
production between the WHIZARD stage and the par- 
ton shower stage. This procedure treats multi-gluon ra- 
diation only approximately and ignores quantum inter- 
ference between the electroweak and QCD production 
amplitudes. However, these are relatively small effects 
at the ILC and are unimportant except in dedicated QCD 
studies. 

About 7 million events were created and processed 
through the full GEANT 4"'^ detector simulation, with in- 
dividual events weighted to reflect the statistical sam- 
pling. However all of the 6 and 8 fermion states, the ones 
most relevant for the analysis, were left unweighted. The 



sample has been subsequently divided into bbqqqq final 
states, which constituted the signal and all remaining 
events representing the background. 

In addition to this 'pseudo data' events a further inde- 
pendent sample of 2 million bbffff events was used for 
the calibration of algorithms. 

The jet clustering algorithm used in this analysis is 
the y-cut algorithm^® with the number of jets fixed at six 
to match the number of jets expected for a hadronic ti 
event. 



V. TOP QUARK SELECTION AND e+e" -^ tt ^ bbqqqq 
CROSS SECTION 

The analysis starts with a simple event selection 
based on several global variables described below. 
Events with isolated leptons, defined as a jet containing 
only one reconstructed particle which iseither an elec- 
tron or a muon are rejected as only the bbqqqq i\na\ state 
was considered. 

Subsequently a set of kinematic and topological dis- 
criminating variables has been defined: the total energy 
of the event; the jet finder y^e parameter, which repre- 
sents the y-cut separation between the five and six jet 
hypothesis; the number of particles and the number of 
tracks. The number of particles in the event is defined 
as the number of reconstructed particles identified by 
the PFA. Figure[2]shows distributions of these variables 
for the signal and background samples before any se- 
lections. 





FIG. 2: Kinematic and Topological Event Selections a) yse, b) 
total energy, c) number of particles in the event, d) number of 
tracks in the event. 



Table [D presents the kinematic and topological event 
selections that have been used. After this stage all but 
492000 background events have been rejected. This 
compares to the initial number of 12.5 x 10^ events. 
The efficiency loss for the initial 143000 signal events 
due to this procedure is equal to 9.7%. The subsequent 



Etot 


> 


400 GeV 


log(y56) 


> 


-8.5 


number of particles in event 


> 


80 


number of tracks in event 


> 


30 



TABLE I: List of the kinematic and topological event selections. 



stage of the analysis aims to identify the b quarks and to 
identify the W bosons exploiting its significant invariant 
mass. 

For the purpose of b quark identification the output of 
LCFI flavour tagging algorithm^^ has been used with the 
default settings. Figure |3]shows the performance of the 
LCFI b-tagging algorithm when used for a e+e^ -^ ti ^ 
bbqqqq sample. The neural network output for uds, c and 
b quark jets demonstrates a good separation of differ- 
ent quark flavours for multi-jet environment. In numerical 
terms, a selection corresponding to the b quark tagging 
efficiency of 45.0% will tag 2.6% of charm quarks and 
0.8% of light quarks^. 
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FIG. 3: Distribution of the flavour tagging neural network output 
for uds, c and b quark jets. 



Three additional event selections have been applied 
to the remaining events. The sum of the b-tag neural 
network outputs of all six jets has been required to be 
higher than 1.5; the b-tag parameter of the most b-like 
and second most b-like jet has been required to be at 
least 0.9 and 0.4 respectively. Figure |4] shows the sum 
of the b-tag of the neural network outputs of all six jets 
for the signal and background events after the kinematic 
and topological event selection and before any b-tagging 
selection. It is clear that this is a powerful discriminant 
to select a clean tt sample. In order to identify the in- 
variant mass of the reconstructed W bosons the KinFit 
kinematic fitting algorithm^^ has been used with a single 
constraint that the masses of the two W bosons were 
equal. Only the four least b-like jets have been consid- 
ered for the fit in order to reduce the number of combi- 
nations. All the events with a W mass of more than 1 1 



Hadronic ttbar 
SM Background 




the e+e -^ tt ^ bbqqqq process by the simple use of 
the following formula: 



4 5 6 

Sum b-tag NN output 



FIG. 4: The sum of the b-tag neural network outputs of all six 
jets for the signal and the background events after the kine- 
matic and topological event selection. 



GeV or less than 50 GeV have been rejected. 

After the b-quark and W boson identification pro- 
cedure approximately 74000 bbqqqq signal and 33500 
background events have passed all selections corre- 
sponding to signal efficiency of 51 .5% and purity of 
68.8%. A significant proportion of the remaining back- 
ground derives from the W+W~ -j> qqqq and bbli^qq 
events, with a smaller contribution from ZZ -^ qqqq. 

The top quark mass has been reconstructed using the 
same kinematic fitting approach. The primary aim of this 
procedure was to find a correct match of the b quarks 
to the corresponding W boson, which will be required 
when the polar angle of the top quark needs to be re- 
constructed later in the analysis. The reconstructed top 
mass was also used to further suppress the background 
rejecting all events with masses lower than 145 GeV and 
higher than 195 GeV. Events that yield a probability of 
less than 1% with respect to the constrains used for the 
fitting are also rejected. All constraints used for the top 
mass kinematic fitting can be found in Table [IB The final 
efficiency of the wholes selection process is 29.8% for a 
purity of 79.7%. 



Na 



LL 



Nbg 



■jLdt 



(2) 



where Nall is the total number of observed events, 
while Nbg is the number of simulated background 
events, e is the signal selection efficiency and / Ldt is 
the integrated luminosity.Under the assumption that the 
signal efficiency and the integrated luminosity can be de- 
termined with negligible errors and that the background 
can be reliably determined and subtracted the statistical 
error on the cross section is equal to ^/Nall / {(■ j Ldt) . 
The cross section has been calculated to be 287.4 ± 
1 .3 fb for the whole sample, 370.5 ± 1 .6 fb and 204.3 ± 
1 .2 fb for the two different polarization samples; the first 
cross section being for the sample with negative elec- 
tron polarization. It has to be noticed that these are the 
cross sections for the e+e" -j> it -j> bbqqqq process and 
not for tt production. 



VI. QUARK CHARGE RECONSTRUCTION 

Next step in the analysis is reconstruction of the quark 
charge which is necessary to determine the forward 
backward asymmetry of the bottom and top quarks. 
Hadronization and fragmentation processes obscure the 
quark charge since the bottom quarks fragment into 
neutral mesons in more than 50% of the cases. While 
charged B mesons, when reconstructed correctly, allow 
for unambiguous interpretation of the quark charge, for 
the neutral B hadrons the charge is not representative 
of the quark charge. Moreover the neutral B mesons 
oscillate which further dilutes the charge reconstruction. 

Several variables sensitive to the charge have been 
studied and an efficient quark charge estimator has 
been devised as a combination of two variables, the ver- 
tex charge and jet charge, as described below. Note 
that this technique considerably improves a simple ver- 
tex charge algorithm used in the LGFI Vertex software^. 



Mass(topl) 


= 


Mass(top2) 


Mass{W1) 


= 


80.4 GeV 


Mass{W2) 


= 


80.4 GeV 


Mass{b.jeti) 


= 


5.8 GeV 


Mass{bjet2) 


= 


5.8 GeV 


Erotai 


= 


500 GeV 


p^;py;pz 


= 






TABLE II: List of kinematic fitting constraints used for the cal- 
culation of the top mass. 



Once the event selection has been performed it is rel- 
atively straightforward to calculate the cross section of 



A. Vertex Charge and Jet Charge Algorithms 

The vertex charge algorithm uses all tracks associ- 
ated to a secondary vertex weighted by their momentum 
to define the vertex charge Qvtx as per the following 
formula: 



iVTX 






A 



(3) 



where Qj is the charge of the j-th track, -pj is the mo- 
mentum of the track and fc is a user defined parameter; 
the sums are performed only on the tracks associated 
with the vertex. The fc parameter was chosen at 0.3 after 



optimization. Tine performance of sucin metinod for dis- 
criminating tine pa rton c liarge in tine signal sample can 
be seen in Figure [5(a)| Only genuine, identified at the 
MC level b quark jets with a neural net b-tag higher than 
0.4 were included without any requirement to the b quark 
final state. 

Another method of the quark charge determination 
implemented in the analysis, momentum weighted jet 
charge^^, is similar to the one already described in 
Equation |3]with the only difference in the track selec- 
tion process which now includes all the tracks present 
in a jet rather than in a vertex. The jet charge algorithm 
recovers 3.2% of identified b-jets which do not have a 
secondary vertex. 
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FIG. 5: Distributions of reconstructed charge for the template 
signal sample for b quark and b quark jets a) using the momen- 
tum weighted vertex charge, b) using the momentum weighted 
jet charge. 



T he pe rformance of the algorithm can be seen in Fig- 
5(b)| Also in this case the optimal value for the k 



ure 

parameter has been determined to be 0.3. 

The two algorithms rely on different principles to iden- 
tify the quark charge. The jet charge algorithm ex- 
ploits the kinematic consideration that the most ener- 
getic hadrons have a higher probability of containing the 
charge of the quark that initiated the jet^. The principle 
behind the vertex charge algorithm is instead based on 



precisely determining all the tracks that derive from the 
displaced vertex due to the b quark considerable lifetime. 
In this case the aim is to directly determine the charge 
of the meson while the momentum weighting folds in in- 
formation on the reliability of the track. 
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FIG. 6: Performance of a) momentum weighted vertex charge 
and b) momentum weighted jet charge in distinguishing B+ 
from B', BO from j5°, and b from b. 



It is expected that the vertex charge algorithm is 
sensitive only to the charged B-mesons while the jet 
charge algorithm is more universal for different B-hadron 
species. The performance of the algorithms has there- 
fore been tested for only charged mesons B+ and B^ 
and for only neutral mesons B" and 5°. Figure |6(a)| 
shows the purity of a sample with a certain quark charge 
as a function of selection efficiency. This demonstrates 
that the momentum weighted vertex charge is able to 
distinguish well between B+ and B^, while having al- 
niost no discriminatory power when it comes to B° and 
BO. 

Differently the performa nce of the momentum 

is more similar be- 



weighted jet charge. Figure 6(b) 



tween the two cases and the algorithm can separate 
reasonably well also B° and B°, which included both 
B^ and 5° mesons. Most of the difference between 
the charged and neutral mesons in this case can be 
attributed to the flavour oscillations of neutral mesons. 



While in tinis process tine cinarge of tine meson does not 
cinange, tine cinarge of tine b quark does. Tinis introduces 
a furtlner dilution in the charge discrimination. The effect 
is rather small in the B'^ mesons, which have a period 
of oscillation larger than their mean lifetime. In the case 
of B° mesons the effect is dominant as oscillations are 
much faster than the meson lifetime. 



B. Combined Charge 

As the two different methods rely on different informa- 
tion and are rather independent, they have been com- 
bined into a single discriminant, based on the probabil- 
ity ratios^^. If fi{xi) is the probability density function for 
the b quark for variable x, and .ff{xi) is the equivalent 
distribution for the b quark then for each discriminating 
variable x, their ratio r, is defined as: 



riixi 






(4) 



where the index i denotes the discriminating variable. 
Distributions of /'' and /^ were determined using inde- 
pendent samples. 

For each data event a combined tagging variable can 
then be defined: 



n^ 



(5) 



The range of possible values for r is between and cx). 
Given the definition of r, if r < 1 then the reconstructed 
jet is more likely to be from a b quark and if r > 1 the 
jet is more likely to originate from b quark. For conve- 
nience a variable C changing between -1 and -i-1 has 
been defined as: 
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FIG. 7: Combined charge a) distributions for 6 quark and h 
quark jets b) purity versus efficiency curves for h quark and h 
quark jets for combined charge, momentum weighted vertex 
charge and momentum weighted jet charge. Shown for the 
1 74.0 GeV sample after all event selections have been applied. 



VII. QUARK FORWARD BACKWARD ASYMMETRIES 



A. Bottom Quark Forward Backward Asymmetry 



C = 



1 



1 



(6) 



A jet with C > is more likely to derive from a b quark 
and a jet with C < is more likely to derive from a b 
quark. Figure [7(a)] shows the combined quark charge 
performance for the 174.0 GeV top quark sam ple af ter 
all event selections have been applied. Figure [7(b) 



in- 



stead shows the purity versus efficiency curves for the 
combined charge algorithm in the same sample when 
compared to standalone momentum weighted vertex 
charge and momentum weighted jet charge algorithms, 
as it can be seen in fig. [7(b)[ The algorithm efficiency is 



improved by 4% to 10% for a purity range from 60% to 
80%. 

The method described above allows a straightforward 
inclusion of other quark charge discriminants such as 
the lepton charge^ and dipole charge^. 



Before calculating the forward backward asymmetries 
for the b and the t quark, the possibility of performing 
an event selection based on the reconstructed charge of 
the quarks has been investigated. For this purpose one 
would like to use the information derived from both jets. 
Assuming that the event is actually a bbqqqq, rather than 
an event from the SM background, and that the quark 
identification has been correctly performed, the charge 
calculations performed on the two b-jets are really two 
uncorrelated measurements of the same quantity. The 
two b-jets must, in fact, have opposite absolute values 
for their charge. 

The combined charges of the two jets with the highest 
neural net b-tags are therefore multiplied and used as 
an event selection parameter. Figure [8]shows such dis- 
tribution for the signal events with explicit contributions 
from mis-identified events where the mis-identification 
occurred either in the b-tagging (mistagged) or in the 



quark charge determination (wrong cinarge). Tine main 
aim of tinis procedure would be not to reject tine SM back- 
ground but ratiner to suppress tine events in winicin tine b 
quark lias been mistagged or tine cinarge of sucin quark 
lias been misreconstructed. An event cinarge is labeled 
as misreconstructed when the reconstructed combined 
charge of the b quark jet is higher than the combined 
charge of the b quark jet. 



exceeds 60%, while in the 9 = -1 region it is only 15%. 
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FIG. 8: (Combined charge b-jet1) x (Combined charge b-jet2) 
distribution for reconstructed events, mistagged evens and 
events with misidentified charge. Only bbqqqq signal events 
are used. 



A n optim ization has been attempted and the value of 
S/yfS + B has been maximized where S is the num- 
ber of signal events and B is the number of back- 
ground events including both the SM background and 
mistagged events. Interestingly enough the optimiza- 
tion suggested that all events should be included. Un- 
der these conditions the total signal efficiency is 22.7%, 
while the signal purity is 58.1%. The impurities derive 
45.9% from the SM background, 45.0% from the charge 
misreconstruction and 9.1% from the b quark misidenti- 
fication. 

The calculation of the forward backward asymmetry 
as defined in Equation [T| can now be performed using 
two jets with the highest values of neural net b-tags. The 
jet with a higher combined charge has been declared as 
originating from a b quark, while the other b jet has been 
declared as originating from a b quark. The angle 9 of 
the reconstructed &jet has been used as an approxima- 
tion the original b quark angle. Figure[9]shows the event 
distribution with respect to cos(6') of the signal and back- 
ground events after all selections. The mistagged and 
SM backgrounds peak in the forward regions where the 
asymmetry is maximal. This emphasizes importance of 
the forward region in the detector design considerations. 
Note that the mistagged events in the distribution include 
a contribution from b quarks which peaks a\6 = -I. This 
explains relatively high mistagging rate at 6* = -1 when 
compared to the number of b quarks from the signal. It is 
because of this reason that in the 6* = 1 region the purity 
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FIG. 9: Number of events used for the calculation of the b quark 
Afb as function of the b quark 6 angle. In order to qualify as a 
b rather than b quark, the combined charge of the jet must be 
higher than the one of the other b jet present in the event. The 
mistagging refers to both quark charge and flavour. 



The Afb calculation proceeds as follows. The num- 
ber of correctly reconstructed bbqqqq events is evalu- 
ated for the forward and backward hemispheres inde- 
pendently. For this the SM background is subtracted 
from the total number of reconstructed events. The num- 
ber of events left is then multiplied by the purity of the 
reconstruction, accounting for all the events where the 
charge has been misidentified or where the b-jet has 
been mistagged. The number of correctly identified b- 
jets is: iVf, = {Ntot - Nsm) * p, where Ntot is the total 
number of reconstructed events, Nsm is the SM back- 
ground and p is the purity of the reconstruction. This 
equation is applied to each hemisphere, correspond- 
ingly separate purities have been calculated for the for- 
ward and the backward hemispheres. In principle, the 
number of events should also be corrected for the signal 
efficiency because Equation []] uses the cross sections. 
However the efficiencies in the forward and backward 
regions to a good approximation cancel each other out 
and produce a negligible effect on the final result.This 
also leads to robustness of the measurement to varia- 
tions of fragmentation and hadronization models which 
results in a negligible systematic uncertainty. Table [Ml] 
shows the Aft results for different event selections. The 
first line corresponds to the case of no selection which 
maximizes the sensitivity as discussed above^''"^°. 

For each event selection the uncertainty has been 
calculated with three different assumptions. The low- 
est uncertainty, cti, assumes that the efficiency of tag- 
ging and the standard model background have been per- 
fectly simulated at the MC level and therefore do not 
contribute to the uncertainty of the forward backward 
asymmetry. The only uncertainty contribution therefore 



Event Selection 


Aft 


CTl 


0"2 


0-3 


Charge bi x Charge b2 < 1 .0 


0.293 


0.006 


0.007 


0.008 


Charge bi x Charge b-z < 0.5 


0.293 


0.006 


0.007 


0.008 


Charge bi x Charge b-i < 0.0 


0.289 


0.007 


0.008 


0.009 



TABLE III: Reconstructed Aft for the b quark and the respec- 
tive uncertainties. The different uncertainties {ai,o-2,CT3) have 
been calculated with different assumptions as explained in the 
text. 



Event Selection 


Aft 


o"i 


0"2 


0-3 


Charge bi x Charge b2 < 1 .0 


0.356 


0.006 


0.007 


0.008 


Charge bi x Charge b2 < 0.5 


0.348 


0.006 


0.007 


0.008 


Charge bi x Charge b2 < 0.0 


0.353 


0.007 


0.008 


0.009 



TABLE IV: Reconstructed Aft for the t quark and the respec- 
tive uncertainties. The different uncertainties (ai, 0-2,0-3) have 
been calculated with different assumptions as explained in the 
text. 



is yiv. 



tot, e<(>)90= 



where K 



tot, e<(>)90' 



ber of events with b quarks reconstructed in the forward 
(backward) region of the detector. For the second eval- 
uation, (72, the statistical uncertainty from the b-tagging 
efficiency is added in quadrature to the value of cti. The 
added statistical uncertainty is calculated from the pre- 
viously mentioned ad-hock generaded calibration sam- 
ple. Finally the third evaluation, 0-3, considers also an 
additional contribution from the statistics of background 
samplesil, which is added in quadrature to 0-2- In each 
of the three cases the uncertainty has been calculated 
separately for the forward and backward regions and 
subsequently the standard error propagation has been 
used to evaluate the Aft uncertainty. 

The calculated asymmetry agrees well with initial 
asymmetry at the MC level, 0.291 , which suggests that 
the performed analysis has not introduced any system- 
atic bias. In order to check for any significant detector 
smearing leading to systematic effects in Afb the angu- 
lar resolution of the b jet 6 angle with respect to the orig- 
inal b quark has been determined. The resolution has 
been found to be 0.08 radians, and therefore its effect 
on the reconstructed asymmetry is negligible. 

Finally the result can also be decomposed with re- 
spect to the different polarizations used. In the case of 
-80% electron polarization and -i-30% positron polariza- 
tion the asymmetry has been calculated to be 0.356 with 
an uncertainty of 0.010. In the case -i-80% electron po- 
larization and -30% positron polarization the asymmetry 
has been calculated to be 0.155 with an uncertainty of 
0.012. In both cases the 0-3 definition of error is being 
used. 



B. Top Quark Forward Backward Asymmetry 

The analysis of the top quark asymmetry is similar 
to the one already presented for the b quark. The only 
added complication is that, differently from the b quark, 
where the angle 9 of the 5 jet can be used as a very good 
approximation to the angle 9 of the original b quark, the 
direction of the top quark must be reconstructed from its 
decay products, using the kinematic fitter to determine 
correct pairing of two b quarks and two W bosons. More 
specifically the direction of the top quark is calculated 
from the combination of jets that minimizes the x^ of the 



is the total num- fit given the constraints stated in Table EB 
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FIG. 10: Number of events used for the calculation of the top 
quark Aft as function of the top quark 9 angle. In order to 
qualify as a i rather than i quark, the combined charge of the 
b-jet used to reconstruct the top quark must be higher than the 
charge of the other b-]et present in the event. 



The charge of the top quark is determined through 
the charge of the daughter b quark. If a reconstructed 
b quark jet is part of the three jets used to reconstruct 
the top quark then the top quark is declared as a t. If 
instead a b jet is present the quark is declared as a t. 
Given the constrains set to the kinematic fitter only one 
such quark will be present in each jet. Figure [To] shows 
the distribution of top quark events with respect to their 
cos(6'). The distribution includes the SM and mistagging 
backgrounds. 

Subsequently the same Aft, calculations have been 
performed as the ones described in the previous section 
for the b quark case with results shown in Table [TVl The 
calculated asymmetry agrees well with the initial one at 
the MC level, 0.351. 

Finally, in the same fashion as for the b quark, the 9 
angle resolution has been found equal to 0.19 radians. 
This will have a negligible contribution to the total cal- 
culated Afb because only the very central events of the 
Figure [To] distribution will ever be smeared enough to 
change hemisphere when reconstructed. The asymme- 



try in this region is inowever small and does not affect the 
total Afb. 

Similary to the asymmetry of the bottom quark the 
achievable statistical precision has been calculated also 
for the two cases of polarized beams. Unsurprisingly 
values identical to the ones presented for the bottom 
quark (0.010 and 0.012) have been found. 



VIII. DISCUSSION 

To put these results in context, we will interpret them 
in terms of constraints on the couplings of the top quark 
to the vector bosons, ttz and Wtb. 

As we have pointed out already in the Introduc- 
tion, many models of new physics predict large correc- 
tions to the left- and right-handed vector ttZ couplings. 
The measurement of forward-backward asymmetries in 
e+e^ -^ ti will allow these couplings to be determined 
experimentally in a very direct way. 

As a starting point for the analysis, we define the 7 
and Z vertex form factors of the top quark by 



*-' + + 



3A^tl7^(PL-FL7 + PRFRn) + i 



2m* 



F2^]t 



+eZ^J[r{PLFLZ + PrFbz) + i 



2mt 



F2ZW) 



where Pl and Pr are the left and right handed chiral 
projection operators, q,, is the 4-momentum of the virtual 
photon or Z° and tr^'" = i/2{-^^'Y-l''l^)■ The tree-level 
Standard Model values of the form factors are 



Ft,-, — Fr^ — - 



-L7 



i?.7 



-F27 — F2Z — 



F, 



LZ 



V2 -j wJ 



F, 



RZ 



(-44) 



(8) 



where s^ and c„, are the sine and cosine of the Wein- 
berg angle. For later reference, the numerical values 
of the Standard Model Z boson form factors, using 



0.231, are 



F, 



LZ 



0.821 



F, 



RZ 



-0.365 



In principle, we could also introduce in each line a 
fourth, CP-violating, form factor proportional to (7^^''y\ 
One might also include contact interactions between the 
e+er and tf states^V 

In principle, a complete helicity analysis of the full set 
of production and decay angles has the power to con- 
strain many of these parameters independently. How- 
ever, in this paper, we have concentrated on the exper- 
imental measurement of the forward-backward asym- 
metries. Since our main concern here is to illus- 
trate the power of that measurement, we will choose a 
parametrization with two free parameters that can be de- 
termined in terms of the two top quark forward-backward 
asymmetries corresponding to two cases of beam polar- 
izations. 



In the following, then, we will assume that the ^ti 
form factors take their Standard Model values given in 
(|8), that the magnetic moment Z form factor F2Z is 
zero, and that the decay form factors take their Stan- 
dard Model values. We will allow only values of the Zti 
form factors Flz and Fb,z to be varied, and we will deter- 
mine these parameters from two measurements of the ti 
forward-backward asymmetry with different beam con- 
ditions. The choice of -80% electron polarization and 
H-30% positron polarization leads to tt production domi- 
nantly from the initial state e^e^. In the Standard Model, 
for this polarization choice, the 7 and Z s-channel ampli- 
tudes interfere constructively for the production of tLtR 
and destructively for the production of tntL, leading to a 
large positive forward-backward asymmetry. The main 
effect of changing the Z form factors is to relax the de- 
structive interference in the production of tfttL- Thus, 
the asymmetry in this polarization state is mainly sen- 
sitive to Fez, which gives the larger effect on the <_r<l 
state. Similarly, the choice of -i-80% electron polarization 
and -30% positron polarization leads to tt production 
dominantly from the initial state e^ej. In the Standard 
Model, for this polarization choice, the 7 and Z s-channel 
amplitudes interfere constructively for the production of 
t^tL and destructively for the production of tLtn. This 
also leads to a large positive forward-backward asym- 
metry, but one that is mainly sensitive to Flz- Thus, 
the measurement of the tt forward-backward asymme- 
try with these two beam settings sensitively picks out 
non-Standard contributions to the two separate Zif vec- 
tor form facto rs^^. 

For 100% polarized beams, the sensitivity of the tt 
forward-backward asymmetries to deviations of the Z 
form factors from their Standard Model values can be 
computed to be 



5Afb{LR)\ 
5Afb{RL) 



0.138 -0.392 
0.461 -0.106 



SFlz 
SFrz 



(9) 



using Vs = 500 GeV and sin^0„, = 0.231. The large 
off-diagonal terms in this matrix show clearly the effect 
discussed in the previous paragraph. For an electron po- 
larization of -80% and a positron polarization of -i-30%, 
the fraction of events in the two relevant initial polariza- 
tion states is 

/,.-4,.<i^£(£:M^i!fe;»= 0.035 ,10) 

By taking this into account, it is possible to transform 
the matrix presented in Equation [1] in order to account 
for the beam polarizations actually used. Recomputing 
the numerator and denominator for Afb, we find that the 
relation between the form factor deviations becomes 



fSAFBiLR)\ 
{5Afb{RL)) 



0.164 
0.367 



:S)(S:)- <"' 
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Then the standard uncertainties reported in Section VII, 

criAFBiLR)) = 0.010 (7{AFBiRL)) = 0.012 (12) 

gives the uncertainties on SFlz and SFjiz, 

cr{SFLz) = 0.051 cr{6Fnz) = 0.042 (13) 

with some correlation between the values. Normalizing 
to the Standard Model values of these parameters, 

a{SFLz)/FLz ^0.062 aiSFnz) / Frz ^ 0-116 (14) 

These uncertainties are comparable to the values sug- 
gested in2^ on the basis of parametric simulations. One 
can see, for example, by comparing the models dis- 
cussed in^°, that such measurements would cut deeply 
into the space of deviations predicted in models of new 
physics. 

The ILC study of the reaction e+e^ -^ tt will also in- 
clude events in which either the t or the i decays lep- 
tonically. These events add a data set of approximately 
equal size to the one considered here in which the t/t 
charge discrimination is unambiguous. Thus, the full 
analysis of the ILC data will do even better at determin- 
ing the Zit couplings. 

In a similar manner the results can also be interpreted 
with respect to the Wtb anomalous couplings. As a mat- 
ter of fact, the decay form factors of the top quark are 
already constrained at the 20% level by the measure- 
ment of the W helicity at hadron colliders- , and these 
measurements will be improved at the LHC. Thus, it is 
likely that, by the time the ILC operates, the decay form 
factors could be fixed to experimentally determined val- 
ues. Nevertheless, for completeness, we consider the 
effects of these anomalous couplings, following the no- 
tation in'*. 

In this case the appropriate vertex under considera- 
tion is: 



-w 



9 



,, - -^[W-b [i^^AlPl + I.ArPb) t 

-Wf,M^''{BLPR. + BRPL)t] (15) 



2M, 



w 



D, 



= d^ - ieAf_,. 



A, 



where VF^^ = D^^W^ - D^Wf, , i^^ = c^ - ze/i^. /±l,r. 
and Bl^r are the coupling form factors. In the Standard 
Model Al is equal to one, while all the other form factors 
are equal to zero. 

Table [ylpresents predictions of the b quark asymmetry 
for different values of the Wtb anomalous couplings'*. 

It can be inferred from Table [V] that the measure- 
ment of the the b quark forward-backward asymmetry is 
sensitive to the presence of a Bl anomalous form fac- 
tor whose absolute value is greater than approximately 
0.05. Measurement of other observables, not consid- 
ered in this paper, that specifically target the top quark 
decay properties will put much stronger constraints on 
both Bl and Br. 

Note that there is a difference between the asymmetry 
predicted by the Whizard generator which was used for 



Br 


Bl 


Afb 


0.0 


0.0 


0.279 


0.0 


-0.2 


0.243 


0.0 


-0.4 


0.218 


0.0 


-0.6 


0.197 


0.0 


-1.0 


0.169 


-0.6 


0.0 


0.301 


-1.0 


0.0 


0.315 



TABLE V: Aft asymmetry of 6 quark from the top decay for the 
Standard Model and anomalous Wtb vertices. Calculated at a 
centre of mass energy of 500 GeV and in the centre of mass 
rest frame. 



these studies, 0.291 , and the asymmetry by the Com- 
pHEP MC generator^'*?^ used when calculating the the- 
oretical predictions'*, 0.279. Part of the discrepancy can 
also be explained by the fact that the generated sig- 
nal sample is an all inclusive e+e~ -^ bqqqq rather than 
e+e" -^ tt ^ bbqqqq as assumed in the theoretical pa- 
per. In any event, this difference is not significant for the 
purpose of sensitivity estimation. 



IX. CONCLUSION 

We therefore conclude that the achievable resolution 
for the forward backward asymmetry of the top quark at 
the ILC in the e+e^ -^ tt ^ bbqqqq channel is approx- 
imately 0.008 for a total luminosity of 500 fb^^ Simi- 
larly the achievable resolution for the b quark resulting 
from the top decay is also 0.008. In the case of polar- 
ized beams the achievable resolution for both the top 
and bottom quark asymmetries is 0.010 and 0.012 for 
the -80% electron polarization, -i-30% positron polariza- 
tion and the -i-80% electron polarization, -30% positron 
polarization respectively. This result allows to constrain 
the theoretically predicted deviations from the Standard 
Model in the presence of an anomalous coupling of the 
Zti and Wtb vertices. In the case of the Ztt coupling the 
resolution on the predicted Standard Model form factor 
is of the order of 0.05 and 0.04 for the Flz and the Frz 
couplings respectively. In the case of the Wtb the per- 
formed analysis is sensitive to the presence of an Bl 
anomalous form factor greater that approximately 0.05. 
The analysis employed realistic detector simulations and 
advanced reconstruction algorithms in the framework of 
the Silicon Detector concept. A new quark charge recon- 
struction algorithm used to discriminate between bottom 
quarks and their anti-quarks allowed to achieve a selec- 
tion purity of up to 80% for an efficiency of about 60%. 
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